Universal pair-breaking in transition metal-substituted iron-pnictide superconductors 
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The experimental transport scattering rate was determined for a wide range of optimally doped 
transition metal-substituted FeAs-based compounds with the ThCr2Si2 (122) crystal structure. The 
maximum transition temperature T c for several Ba-, Sr-, and Ca-based 122 systems follows a univer- 
sal rate of suppression with increasing scattering rate indicative of a common pair-breaking mech- 
anism. Extraction of standard pair-breaking parameters puts a limit of ~26 K on the maximum 
T c for all transition metal-substituted 122 systems, in agreement with experimental observations, 
and sets a critical scattering rate of 1.5 x 10 14 s _1 for the suppression of the superconducting phase. 
The observed critical scattering rate is much weaker than that expected for a sign-changing order 
parameter with strong interband scattering, providing important constraints on the nature of the 
superconducting gap in the 122 family of iron-based superconductors. 
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The discovery of iron-based superconductors in 2008 
breathed new life into the study of high temperature su- 
perconductivity, with numerous families of compounds 
since discovered, characterized and extensively studied 
|l|. Intermetallic iron-based systems with the ThCr2Si2 
"122" structure and doped with transition metal (TM) 
substitution on the iron site remain the most widely stud- 
ied due to the feasibility of synthesizing large single-phase 
crystals coupled with the ability to substitute a plethora 
of TM elements for iron. With superconductivity induced 
by substuting almost any of the TM elements in the Fe, 
Co, and Ni columns, the robustness of these superconduc- 
tors to disorder - in particular, disorder focused directly 
in the active pairing layer - provides a striking constrast 
to the sensitivity found in other unconventional super- 
conductors. Furthermore, this robustness initiated one 
of the early challenges to the proposed s± sign-changing 
gap symmetry [2j, and has been toted as evidence for a 
non-sign-changing s-wave pairing symmetry [3|. 

The role of TM substitution in both promoting a su- 
perconducting state and shaping the phase diagrams of 
the 122 systems is an important topic of ongoing debate. 
In the Ba-based 122 systems, the substitution-induced 
positioning of the superconducting phase scale reason- 
ably well with d-electron count (with the exception of Cu 
substitution) [4j] , and ample evidence of modifications to 
band structure, carrier concentrations and magnetic in- 
teractions support a rigid band shift doping model [l|. 
However, theoretical models predicting the localization 
of added <i-electrons and the importance of impurities 
raise questions about this approach Jala]- Moreover, the 
similarity of the phase diagram produced by nominally 
isovalent Ru substitution [7|, |8(] to that of its aliovalent 
counterparts necessitates a better understanding of the 
true nature of TM substitution. 

In this study we compare the elastic transport scat- 
tering rate and superconducting transition temperature 



T c for a wide range of optimally doped TM-substituted 
122 compounds and observe a universal correlation that 
follows an Abrikosov-Gor'kov (AG)-like pair-breaking 
suppression for all types of transition metal substituents 
and alkaline earth cations. We show that the large vari- 
ations in optimal T c values found in different 122 sys- 
tems are due to variations in impurity scattering rate, 
but are also limited by an ideal zero-scattering limit that 
lies much below the T c values of alkali metal-doped 122 
systems. We discuss implications for order parameter 
structure as well as constraints on the inter- and intra- 
band coupling strength determined by the universal re- 
lation. 

Single crystals of AFe2_ a; TM; C As2 compounds (with 
A=Ba, Sr; TM=Co, Ni, Pd, Pt) were synthesized us- 
ing the FeAs self- flux method described previously [9|. 
TM concentrations were determined by wavelength dis- 
persive x-ray spectroscopy (WDS). Resistivity and Hall 
effect data were measured in a commercial cryostat with 
Hall coefficient {Rh) values obtained by antisymmetriz- 
ing field sweeps at constant temperature between -5 T 
and +5 T. To minimize geometric factor errors in de- 
termining scattering rates, six-wire measurements were 
used to simultaneously determine longidudinal (p xx ) and 
transverse (p xy ) resistivities, using gold wire and silver 
paint contacts with typical contact resistances of ~ 1 Q. 
Transition temperatures for samples measured by our 
group were determined by the mid-point of the resistive 
transitions, and are well documented in previous pub- 
lications to coincide with magnetization measurements 
9, 12]. Data obtained from the literature utilized the 
same criteria when possible, and used stated values oth- 
erwise. 

Unlike the BaFe2As2 family, the maximum or "opti- 
mal" transition temperature T c i max \ for different TM- 
doped versions of the SrFe2As2 system exhibits a wide 
variation of values, reaching ~20 K for Co, Rh and Ir 
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FIG. 1: Phase diagrams of the SrFe2-iNi I As2 [9| and 
SrFe2- :c Pt :c As2 Q systems. Antiferromagnetic (triangles) 
and superconducting (diamonds) transition temperatures are 
plotted for Ni- (open symbols) and Pt-doped (closed sym- 
bols) systems. The similar rate of suppression of the magnetic 
phases and the position of the superconducting domes, with 
optimal doping at x ~ 0.16 for both cases, is to be contrasted 
with the considerably different T c ( maa .) values of 9 K and 16 K 
for Ni- and Pt-doped series, respectively. 
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12j, 9 K for Ni [9j and 8 K for Pd 



16 K for Pt 

substitution [10j. As a prime example, we directly com- 
pare the phase diag rams of the SrFe2- a :Ni :E As2 [9] and 
SrFe2- a; Pt a ;As2 12[ systems in Fig. [TJ As shown, the 



antiferromagnetic order transition T/v follows an almost 
identical decline as a function of either Pt or Ni substitu- 
tion, with minimal difference between the two systems. 
The similar positioning of the superconducting dome for 
each system at an optimal concentration of x ~ 0.16 fol- 
lows that expected for the nominally equivalent addition 
of two d-electrons from both Pt and Ni substituents, as 
compared to that of SrFe2- a: Co a ;As2 with only one d- 
electron contribution and a significantly larger optimal 
doping of x ~ 0.24 13j. However, a significant factor 
of two difference is apparent in T c ( max j values, present- 
ing an intriguing contrast in two systems with nominally 
identical phase diagrams. 



With similar modification of unit cell parameters [TJ] , 
identical oxidation states and nearly identical phase di- 
agrams in both substitution series, we consider intrinsic 
variations in pair-breaking scattering rates as the pri- 
mary origin of this contrast. Following previous studies, 
which have shown that electron bands dominate trans- 



port in the TM-doped systems [14j-|l6| and optical con- 
ductivity studies which indicate a single dominant Drude- 
like component 17J, Il8j . we utilize a simple one-band 
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FIG. 2: Six-wire measurements of resistivity (main panel) and 
Hall coefficient (inset) of optimally-doped SrFe2- a; Pd :E As2, 
SrFe2-^Nia;As2, SrFe2-iCoi;As2 and BaFe2-iPta;As2, with 
T c values of 7, 8, 17 and 20 K, respectively. The schematic 
depicts the configuration of the six-wire measurement. 



model [19j, |20| to estimate the elastic (T = 0) transport 
scattering rate, allowing for broad comparisons between 
different systems and data sets. The normal state scat- 
tering rate, T = ^Srjf , where e is the electronic charge 
and 77i* is the effective mass, is determined from trans- 
port measurements by extrapolating the resistivity p xx 
and Hall coefficient Rh to zero temperature using power 
law fits over an extended range of temperatures above 
T c . We use an effective mass of m* = 2m e based on the 
measured values for the electron bands in CaFe2As2 [231 ] . 
SrFe2As2 21 1 and BaFe2As2 [22| from quantum oscilla- 
tion measurements and an optical conductivity measure- 
ment on optimally Co-doped BaFe2As2 [24J . 

Six-wire measurements were used to determine both 
p xx and Rh simultaneously for several optimally-doped 
samples with a range of T c values, thereby significantly 
reducing geometric factor error in calculating F by elim- 
inating sample thickness dependence. Shown in Fig. [21 
four different optimally-doped TM-substituted samples 
(Ni-, Pd-, Co- and Pt-based with T c values of 7 K, 8 K, 
17 K, and 20 K, respectively) exhibit an observable dif- 
ference in absolute resistivity values dominated by a rigid 
shift in the zero-temperature elastic contribution po, as 
evident from the comparable inelastic contributions (i.e., 
slope of p(T)). The resulting contrast in T c values fol- 
lows this trend, with a systematic reduction of T c with 
increasing po. 

We compare the resultant T c (r) values with those cal- 
culated for all optimally doped TM-substituted 122 sam- 
ples with p xx and Rh values available in the literature, 



as shown in Fig. [3] (All data corresponds to systems 
with electron-dominated transport with the exception of 
Ru substitution, for which we utilize the electron com- 
ponent of p xx extracted with a two-band analysis and 
n e from ARPES measurements 37| to obtain a value of T 
that can be compared with the other data.) Remarkably, 
all T c values follow the same trend of suppression with in- 
creasing r, as expected in the AG formalism for a super- 
conductor with increasing levels of pair-breaking impuri- 
ties 20, I25U27I ] , but surprising in light of the variety of 
systems presented. In particular, there is no clear trend 
associated with species of alkaline earth cation or tran- 
sition metal substituent except for an average reduced 
scattering rate for Ba-based systems. This is likely cor- 
related with both the lower substitution concentrations 
required to reach optimal doping as well as the lower 
T/v ordering temperatures in BaFe2As2 as compared to 
both SrFe2As2 and CaFe2As2 . Note that T values for 
BaFe2-a;Co x As2 and BaFe2- :r Ni :E As2are nearly identical 
to those obtained in optical conductivity measurements 
[18| if we assume the same effective mass values, provid- 
ing a confirmation of our analysis. 

The rate of suppression of T c , defined by the critical 
scattering rate r c where T c is completely suppressed, 
is in general dependent on the type of scatterers and 
the order parameter symmetry: according to Ander- 
son's theorem, fully-gapped s-wave superconductors only 
respond to magnetic impurities, while unconventional 
pairing symmetries can be affected by both magnetic 
and nonmagnetic impurities [28J. Assuming predomi- 
nant non-magnetic scattering as evidenced by a para- 
magnetic normal state and no obvious indication of en- 
hanced magnetism due to TM substitution (e.g., absence 
of any enhanced susceptibility) [9|, |l2|, the presence of 
non-magnetic pair-breaking points to a sign- changing or- 
der parameter. However, several substitution 27l.l29j| and 
irradiation 2JJ, I30L l31| studies report a much weaker rate 
of suppression than that expected for a sign-changing or- 
der parameter; calculations for an ideal s± superconduc- 
tor with full gaps on both bands [26j and strong interband 
scattering yield T c (s±) = 1.8xl0 13 s" 1 0, El El] , with 
similar values for the d-wave case [2(j| . Shown in Fig. 02 
a fit to the typical AG functional form [25] yields a value 
r c = 1.5 x 10 14 s^ 1 corresponding to a critical mean free 
path of ~ 1.1 nm (using Fermi velocity Vf = 1.7x 10 5 m/s 
[2lJ), close to the exp ected superconducting coherence 
length £ = 2.8 nm 0. 

However, this value is also an order of magnitude 
weaker than the expected r c (s±), presenting a signifi- 
cant challenge to models considering a fully gapped s± 
pairing symmetry, particularly in the presence of strong 
interband scattering [3J]. But calculations using the T- 
matrix approximation for an s± state emphasize that 
both inter- and intra-band scattering in the unitary limit 
can be decreased with appropriate parameters, resulting 
in a possible four- fold increase of T c 34[ that may offer an 



30 



20 



10 



*9? 



T&r^C^ 







AE[Fe,TM] 2 As 2 


A 


BaFe Ru As 

1.3 0.7 2 


A 


BaFe Pt As 

1.90 0.10 2 


■ 


SrFe Co As 

1,7 03 2 


▲ 


SrFe Pt As 

1 84 16 2 


• 


BaFe Co Cu 

1.90 004 o.c 


♦ 


SrFe Ni As 


k. 


SrFe Pd As 

1843 0.157 



i 



r SJ 



r(xio 14 s" 1 ) 



FIG. 3: Effect of transition metal substitution on T c val- 
ues of a wide variety of 122 superconductors at optimal dop- 
ing concentrations, plotted as a function of the experimental 
transport scattering rate F = epjRura" (see text). Closed 
symbols indicate six- wire measurements (see text), and open 
symbols indicate values obtained from literature data for 
(1) CaFei.92Co .o8As 2 [H, (2) BaFez-.Co^Asa Gl, H|, (3) 
BaFe2-xNi,As 2 [|3, S3, (4) SrEei.wCoo.aeAsa g| and (5) 
SrFei.8 4 Nio.i 6 As2 (po=212 /xSlcm, R H (0 K)=1.4xl(T 9 m 3 /C) 
[5(| . Samples of SrFe2- a; Co a ;As2 and SrFe2- a; Ni a; As2 denoted 
with an "a" are annealed (see text and footnote Q). All 
T c values follow a universal rate of suppression with V well 
described by an Abrikosov-Gor'kov fit (dashed line) that is 
much weaker than expected for a superconductor with s± 
symmetry and interband scattering (dotted line) [20l. l3Cl. l32j| . 



explanation, and may in fact be used to extract the rel- 
ative strength of inter- and intraband scattering in these 
systems. In addition, recent studies that consider the ef- 
fects of disorder on both superconductivity and compet- 
ing states [351 ] provide an alternative explanation for the 
apparent weak pair-breaking effects observed throughout 
the iron-based superconductor family. 

The optimum clean-limit (r=0) transition tempera- 
ture T c q is an important parameter since it is the value 
that should be utilized in considering the intrinsic pair- 
ing strength. Our determination of T c q — 26 K is con- 
sistent with the well-established maximum T c value of 
- 25 K found among all TM-doped 122 systems 0, 
as well as with extrapolated estimations of pressure- 
and doping-optimized systems such as shown in the 
comparison of BaFc2- a; Ru a ;As2 substitution and pres- 
sure dependence [36(. But this observation raises an 
intriguing question about why T c q does not approach 
that found in higher T c intermetallic systems including 
Sri_ ;c K a; Fe2As2, Ba i _ a; K a; Fe2As2 and BaFe 2 As 2 - a; Pa ; 0, 



which have calculated T values in the range shown in 
Fig. [3) In contrast to the typical explanation of a reduced 
level of active-plane disorder as the reason for higher 
T c values in the alkali metal-doped systems, the deter- 
mination of T c0 and its failure to reach ~ 40 K suggests 
a fundamental asymmetry in pairing strength between 
electron- and hole-doped systems that does not arise from 
scattering differences alone (althou gh t he effects of strong 
scattering in the hole bands [14], LL5J, [ljj |37j cannot be dis- 
counted as a factor in the observed asymmetry). 

A universal T C (T) relation suggests a similar pairing 
potential for all TM-doped 122 compounds that is dis- 
rupted by a common scattering mechanism. It is not 
clear why certain TM substitutions induce more scat- 
tering than others, but dramatic variations in seem- 
ingly similar elemental substitutions are not unprece- 
dented. For instance, the BaFe2-zRu :z As2 system re- 
quires ~ 30-40% Ru substitution to obtain optimal dop- 
ing, which is almost four times higher concentration than 
Co substitution but results in a very similar value of 
T c (max) ■ Such a contrast has been argued to arise from 
the aliovalent versus isovalent nature of, respectively, 
Co and Ru substituents, but recent work has put this 
into question. Mossbauer studies of BaFe2- :E Co :E As2 and 
BaFe2-xNia;As2 find no change in d-electron population 
with substitution [38JJ , while an x-ray absorption study 
reveals no change in the Fe valence with Co substitution 
in BaFe2- :E Co a As2 [39]. Furthermore, recent calculations 
suggest that substituted (^-electrons can remain localized 
at the substituent sites [5j , either still resulting in a rigid 
band shift [ff or generating a phase diagram strikingly 
similar to that expected from a rigid band shift [40J. 

Variations in impurity or disorder levels due to details 
of substitution chemistry likely play a key role in ex- 
plaining the variation in Y values observed in the 122 
series of superconductors. This is corroborated by ob- 
servations of enhancements in T c values after annealing 
crystals of both low- and high-r systems, in particu- 
lar BaFe2- :r Co : iAs2 [4l| and SrFe2- :E Ni a ;As2 [42J, respec- 
tively, and confirmed by our study of a SrFe2- a; Ni a ;As2 
crystal with six-wire measurements obtained before and 
after annealing: as shown in Fig. |3l the shift of data 
along the AG curve indicates an inverse relation between 
T c and Y [43]. In the case of BaFe 2 - a: Co :r As2 [4l|, an- 
nealing was shown to enhance T c to a maximum value 
of 25 K, consistent with our determined T c q value. The 
reason why the BaFe2_ x Co :z As2 system is closest to the 
clean limit is not known, however a lack of observable 
disorder in Fe-As bond lengths in BaFe2- x Co x As2 may 
have provided an important insight |44|; it would be in- 
teresting to perform the same study on high-r systems 
to confirm this scenario. 

In conclusion, we have demonstrated the existence of 
a universal pair-breaking relation for a wide range of op- 
timally transition metal-doped 122 systems, suggesting 
a common scattering mechanism and pairing potential 



across the series. The rate of suppression of T c and the 
contrast between the optimum (zero-scattering) T c q value 
of ^26 K and the higher T c values achieved in non- 
transition metal substitution series provides important 
constraints on the pairing symmetry and mechanism in 
the intermetallic iron-based superconductors. 
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